The Mixmaster Universe is Unambiguously Chaotic 
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OS : The Mixmaster or Bianchi IX cosmological model has become one of the archetypal 

settings for studying gravitational dynamics. The past decade has seen a vigourous 
debate about whether or not the Mixmaster's dynamics is chaotic. In this talk we 
review our recent work in uncovering a chaotic invariant set of orbits in the Mixmaster 
phase space, and how we used this discovery to prove the dynamics is chaotic. 
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1 The mixmaster universe 

The mixmaster modefli describes a homogeneous but anisotropic closed cosmology 
of Bianchi type IX. The spatial sections have the geometry of a deformed three- 
sphere, evolving with three independent scale factors, a{t),b{t),c{t). The vacuum 



, field equations read: 
O' (Ina^)" = (62 _c2)2 -a* et eye. {a,b,c) . (1) 

. Here a prime denotes d/dr and dt — abc dr. When tlic right hand sides of the three 

' equations (|l]) are all w 0, ie. when the minisuperspace potential 
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y = a* + 54 + - 2{a'b' + a'c' + b'c') « , (2) 

the universe evolves in an approximate Kasner phase with a ^ t^' , b ^ t^^ , c ^ V'' . 
The Kasner exponents satisfy Yl^i=iPi — Si=iPf = 1- 

The dynamics of the model can be separated into the evolution of atL overall 
scale factor J7, and two variables that measure the departure from isotropp: 

n^Un{abc), /3+ = iln(&c/a2), = -L in(6/c) . (3) 

Since it can be showiH that the evolution of describes a smooth overall expansion 
and contraction of the universe, any chaotic behaviour must be associated with 
motion in the anisotropy plane 

Early work on the mixmaster dynamics focused on discrete mapffl that were ua 
to approximate the continuum dynamics. The maps were shown to be chaotid 
and from this it was inferred that the full dynamics was also chaotic. The picture 
became clouded in the late eighties when nunierical studies of the dynamics found 
that the principle Lyapunov exponent vanishecO, indicating that the system was not 
chaotic. Subsequently, two possible causes for the discrepancy were isolatecfl: (1) 
the full dynamics might be integrable, but the approximate maps fail to preserve 
all the constants of motion; (2) the full dynamics is chaotic but the Lyapunov 
exponents vanish for certain choices of time coordinate. The second possibility 
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raised the important issue of how to describe chaos in a coordinate invariant way. 
We address this question in our other contribution to this volume. In the remainder 
of this talk we describe how coordinate independent methods can be used to settle 
the mixmaster debate in favour of option (2). That is, the mixmaster universe is 
chaotic and the approximate maps do capture the essential features of the dynamics. 

2 Uncovering the chaotic invariant set 

The key to understanding the mixmaster dynamics can be found in the very first pa- 
pers on the subjccoB. These papers identify three "attractors" for the dynamics cor- 
responding to the three states of highly anisotropic expansion a ^ b ^ const., c ^ t, 
c ^ a ^ const., 6 ~ t and b c ^ const., a ^ t. The 5*0(3) symmetry of the spatial 
sections ensures that all three possibilities occur with equal probability. If these at- 
tractors corresponded to final outcomes, it would be a simple matter to decide if the 
dynamics was chaotic. All one would have to do is assign a colour to each outcome, 
then colour-code the space of initial conditions according to the outcomes. If the 
boundaries between the three outcomes are fractal, then the dynamics is chaotic. 




Figure 1: (i) The minisuperspace potential with exits cut in the corners and outcomes marked by 
arrows, (ii) Outcome basins and their fractal boundaries in the anisotropy plane. 

Unfortunately, the procedure is not so simple as the attractors do not cor- 
respond to final outcomes. When viewed in the anisotropy plane, the attractors 
correspond to trajectories that go deep into one of the three corners of the min- 
isuperspace potential (see Fig.(l.i)). After a long time, the trajectory returns to 
the scattering region close to the centre of the triangular potential before again 
scattering toward one of the three attractors. 

We reasoned that if the mixmaster was chaotic, it must correspond to a non- 
compact billiarci]. Moreover, there should be a chaotic invariant set of unstable 
periodic orbits in the centre of the scattering region. In order to uncover this set 
we set a threshold based on the relative rate of expansion of the three scale factors. 
When one axis was collapsing very much faster than the other two we stopped 
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the evolution and assigned an outcome. This is equivalent to cutting pockets in the 
corners of the minisuperspace potential, as shown in Fig.(l.i). In effect wepSimplified 
the analysis by converting the mixmaster into a Hamiltonian exit systerrS. We can 
do this with confidence since a system with exits is always less chaotic than the 
corresponding system without exits. 

With exits in place it is a simple matter to numerically plot the attractor 
basin boundaries. A representative slice through phase space is shown in Fig.(I.ii). 
The basiii. boundaries are clearly fractal, and were found to have an information 
dimensioiia of Di = 1.86 ± 0.02. From this we can conclude in a coordinate inde- 
pendent way that the mixmaster universe is chaotic. Moreover, exactly the same 
procedure can be used to show that the Bianchi VIII cosmological model is chaotic. 

Trajectories belonging to the boundaries of the attractor basins are future 
asymptotic to the unstable periodic orbits that comprise the mixmaster's chaotic 
invariant set. (There are actually no truly periodic orbits in the plane since 

the overall anisotropy increases as a cosmological singularity is approached. The 
orbits are self-similar rather than periodic. However, rescaled coordinates can be 
found where the orbits are periodic) . The complexity of the mixmaster system can 
be measured by introducing a symbolic coding for the periodic orbits. (A similar 
coding was introduced by Rugm for the aperiodic orbits). The most efficient way 
of coding an orbit is shown in Fig.(l.i). The symbol x, y, or z is recorded each 
time the trajectory passes through one of the three triangular regions bounded by 
the corners of the potential. The simplest periodic orbit, first found by MisnetJ, 
corresponds to the repeated sequence xyz. All the periodic orbits can be similarly 
described, and it is easy to show that the number of orbits made from words of 
length k grows as N{k) ~ 2'^. Since the number of orbits grows exponentially with 
length, the collection of periodic orbits has a non-zero topological entropy, in this 
case, Ht = In 2, which implies the dynamics is chaotic. 
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